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In one scenario of baryogenesis, the matter-antimatter asymmetry was generated in the early universe during a 
cold electroweak transition. We model this transition by changing the sign of the effective mass-squared parameter 
of the Higgs field from positive to negative. The resulting 'tachyonic' instability leads to a rapid growth of 
occupation numbers, such that a classical approximation can be made in computing subsequent developments in 
real time. 

We solve the classical equations of motion in the SU(2)-Higgs model under the influence of effective CP- 
violation. The resulting baryon asymmetry follows from the generated Chern-Simons number using the anomaly 
equation. The 'classical' difficulties with lattice implementations of these observables are avoided here because 
the fields are smooth on the lattice scale. 



1. Introduction suffers from short-distance fluctuations. Out of 



A basic ingredient in the theories of baryogen- 
esis is the anomaly in the baryon-current diver- 
gence, 



equilibrium the problems are worse in the quan- 
tum world. 



where q is the topological-charge density in the 
SU(2) gauge helds, which is the divergence of the 
Chern-Simons current, q — 9 M Jq S (the • • • de- 
note the U(l) contribution which is supposed to 
be inessential). In electroweak baryogenesis the 
baryon number B is generated during the elec- 
troweak transition and given by 
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Fortunately, the classical approximation can be 
used in case of large occupation numbers, e.g. at 
sufficiently high temperature. But even in this 
case there are short-distance problems because 
of the slow decrease of occupation numbers at 
large momentum, n p ~ T/p. However, if we 
have good reason to start with suppressed short- 
distance modes, then the problem may show up 
only at late times, because classical equipartition 
is a slow process. For times of interest to the 
problem we can then use a 'naive' lattice tran- 
scription of q. An ideal case is the scenario of 
baryogenesis from a rapid cold electroweak tran- 
sition. 



where Acs = / d 3 xjQg is the Chern-Simons num- 
ber. We assumed the transition to start at t = 



2. Scenario [Tl2] 



The computation of q is a 'classical' problem 
in lattice gauge theory. The 'naive' lattice tran- 
scription 



and set A cs (0) = 0. 
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changes sign, e.g. due to a coupling to the in- 
ffaton field. The negative m^ ff causes a ('tachy- 
onic' or 'spinodal') instability in which the sys- 
tem is far out of equilibrium, and under influ- 
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ence of a CP bias one expects a net change in 
the Chern-Simons number, with a correspond- 
ing baryon asymmetry. Non-linear effects lead 
to rapid approximate thermalization to an effec- 
tive temperature T that is expected (counting 
d.o.f.) to be much lower than the standard finite- 
temperature transition temperature T c , so sub- 
sequent sphaleron transitions are expected to be 
negligible. 

The reason for the change of sign of m 2 s is un- 
certain. Not committing ourselves to any specific 
cause, we use a sudden quench: 

t < 0, 
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t > o. 



This maximal non-equilibrium setup presumably 
also maximizes the asymmetry for given CP vio- 
lation. 



so grows large whereas £~ — > 0. The canonical 
coordinates and momenta grow large and become 
correlated, an ideal case for making a classical 
approximation. 

As soon as the unstable modes have grown large 
enough, but before non-linear interaction effects 
are expected to take over, we could sample the 
above distribution to provide initial conditions for 
subsequent classical evolution. However, we may 
as well sample the distribution at time zero, since 
this leads to exactly the same distribution later on 
in the free-field approximation. As n = n = at 
t = 0, only the 1/2 remain in the denominators, 
which is why we call these 'just a half initial con- 
ditions [2] • Note that only momenta with p < \x 
are initialized. 

4. Simulation 



3. Classical approximation 

The effect of the quench can be studied by first 
neglecting all interactions. Let be a real com- 
ponent of the Higgs field operator. Its spatial 
Fourier transform behaves around time zero as 

' ' ' aipe 1 ^ 1 ) , t < 0, 
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Matching at t — determines a and f) in terms 
of a and fit . The modes p < fi are unstable and 
grow exponentially: uj~ = i\u>l\ — i\/p 2 — p 2 , 



So there is classical behavior for \u~\t 1. In- 
deed, generic quantum correlation functions in 
the initial state |0) (annihilated by the a) can be 
well approximated by a classical distribution 2 , 
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where £p = (ui p </> p ± 7r p ) /^/2cj p , and n and h are 
instantaneous particle numbers. In the unstable 
region, as [it 3> f : 



n p + 1/2 + n p > 1, 
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To study the emerging asymmetry we used the 
SU(2) Higgs model with effective CP violation, 
given by 
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Here k parametrizes the strength of an effec- 
tive CP-violating term of dimension six. ft is 
supposed to represent a leading effect induced 
by physics beyond the standard model, or even 
mimic CP violation in the standard model itself. 

We discretized the action on a (lorentzian) 
space-time lattice, from which the field equations 
follow in the usual way. Initial conditions were 
chosen according to the 'just a half scheme, and 
also using a thermal Bose-Einstein distribution 
with small T — 0.1 m# for comparison. Since the 
initial field configurations are smooth on the lat- 
tice scale, we used 'naive' lattice transcriptions, 
even for tr FF, 

The results to follow where obtained with a 
lattice spacing arriH = 0.35, volume (Litih) 3 = 
(60 x 0.35) 3 = 21 3 , couplings g = 2/3 and A 
such that niH/mw = 1, V2 and 2 (recall m# = 
\/2Av = \/2jU, raw = gv/2, with v = p/VX the 



3 



vacuum expectation value of the Higgs field) . 

Figure ^ shows the behavior in time of the 
volume-averaged 0^0, the Chern-Simons num- 
ber and the winding number in the Higgs field, 
A w , for one initial condition, for k — 0. At 
first <jy <fi increases exponentially, then peaks at 
mnt = 7.5 and subsequently executes a damped 
oscillation towards a value somewhat smaller than 
v 2 . For small energies, in equilibrium, the wind- 
ing number is expected to be near Acs, as this 
diminishes the covariant derivative D^(j>. Fig. ^ 
shows that A w behaves at first erratically and 
then settles near an integer (the smaller the lat- 
tice spacing, the closer to an integer). The early 
Acs is very small and subsequently it appears 
to approach A w . We see no sphaleron transi- 
tions (jumps in Acs by approximately an inte- 
ger) in 50 < mnt < 500. The effect of non-zero k 
on single trajectories is unpredictable, due to the 
chaotic nature of the dynamics, but an effect is 
present in the average over initial conditions (as 
needed for the quantum expectation value). 

Figure |21 shows an example of averages for non- 
zero k. After an initial rise of (Acs) (which we 
understand semi analytically), there appears to 
be a sort of resonance with {4>'<p) that may even 
lead (for larger k) to a (Acs) 01 opposite sign. 
We found a similar resonant effect in 1+1 dimen- 
sions, where it leads to a strong dependence on 
mn/mw |2| Note that after tmn — 40 there 
is no visible drift of (Acs) towards zero due to 
sphaleron transitions. 
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Figure 1. Example of cf^(f)/v 2 , A w and Acs versus 
time; mn/'mw = 1, thermal initial condition. 
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Figure 2. (ft <j)/v 2 ) , (Acs) and (A w ) for k = 
(l/27r 2 )m^, 2 , mn/rnw = and 'just a half 
initial conditions. 



5. Conclusion 

Baryogenesis through tachyonic electroweak 
transition offers ideal case for lattice simulations 
using the classical approximation. Together with 
Jon-Ivar Skullerud we have also computed Higgs 
and W particle numbers [jj]. As expected, they 
grow large at low momenta, while staying negligi- 
ble for momenta p > 1.5 Wjj. This supports our 
experience that the effective CP-violating interac- 
tion — Ccp = K J <jr (pti FF can be meaningfully 
implemented using a 'naive' lattice transcription. 
More details are in 0] (and to be published). 

The resulting asymmetry appears to be sub- 
stantial: using simple estimates the observed 
baryon to photon ratio can be reproduced with 
k « 1 x 10~ 5 TeV -2 , assuming m# = \[2m-w- 
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